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-^ . Abstract 

pi i' In the recent years there has been an enormous development in the evaluation of 

^ • higher order quantum corrections. An essential ingredient in the practical calcula- 

tions is provided by vacuum diagrams, i.e. integrals without external momenta. In 
this paper a program package is described which can deal with one-, two- and three- 
^ ■ loop vacuum integrals with one non-zero mass parameter. The principle structure 

is introduced and the main parts of the package are described in detail. Explicit 
examples demonstrate the fields of application. 



PROGRAM SUMMARY 

Title of program: MATAD 

Available from: 

pittp : //www-ttp . physik . uni-karlsruhe . de/Progdata/MATAD/l/| 

Computer for which the program is designed and others on which it is operable: Any 
work-station or PC where FORM is running. 

Operating system or monitor under which the program has been tested: UNIX, 
FORM 2.3 

No. of bytes in distributed program including test data etc.: 706000 

Distribution format: ASCII 

Keywords: three-loop computations, vacuum integrals, computer algebra, automa- 
tion of computations 

Nature of physical problem: Multi-loop integrals are needed for the evaluation of 
quantum corrections. An important class of loop diagrams is covered by so-called 
vacuum integrals which have no external momentum. MATAD can analytically com- 
pute those one-, two- and three- loop vacuum integrals where one mass scale is 
present. 

Method of solution: The method of integration-by-parts is used in order to obtain 
recurrence relations which reduce complicated integrals to a small set of so-called 
master integrals. They have to be evaluated once and for all. In addition a user 
interface is provided which makes it easy to put in complicated diagrams in a rather 
compact way. 

Restrictions on the complexity of the problem: The restrictions on the complexity 
are given by the hardware limitations of the computer and the limits on the size of 
the storage files inside FORM. 

Typical running time: The runtime strongly depends on the complexity of the di- 
agram under consideration. It may vary form a few seconds to the order of a few 

weeks. 



LONG WRITE-UP 



1 Introduction 

The high experimental precision reached at the electron-positron machines LEP (CERN) 
and SLC (SLAC) and the hadron collider TEVATRON (FERMILAB) requires from the 
theoretical side the evaluation of higher order quantum corrections. In the cases where 
perturbative methods are applied the quantum corrections can be expressed through 
an expansion in the coupling constant of the underlying theory. The individual terms 
can in turn be expressed through so-called Feynman diagrams, which are often classified 
as multi-leg or multi-loop diagrams. Vacuum integrals, i.e. integrals without external 
momenta, constitute an important sub-class and often serve as building blocks in complex 
calculations. 

In general the momentum integration of the loop integrals is divergent in four space- 
time dimensions. At present the most practical method to cope with this problem in 
higher loop orders is based on Dimensional Regularization [jl|. There, the four space-time 
dimensions are replaced hj D = 4 — 2e dimensions. Then the integrals are solved for a 
choice of e that renders them finite. Finally an expansion for e — ;► is performed and the 
divergences manifest themselves as poles in e. 

Important progress in practical computations has been made roughly 20 years ago 
by establishing an algorithm for the evaluation of propagator-type diagrams up to three 
loops in the massless case 0] . They are important if there is only one external momentum 
which sets the mass scale for the problem. The formulae have been implemented on 
the computer in the FORM |^ package MINCER [^. In 1995 for the first time three-loop 
diagrams in the opposite limit, i.e. zero external momentum but massive lines, were 
systematically examined [^. Usually these are denoted as vacuum or tadpole diagrams. 
In all integrals contributing to the photon propagator have been considered. The main 
characteristics of this class of diagrams is that the massive line forms a closed loop. These 
considerations have been extended to the W boson current correlators which led to one of 
the most prominent applications of three-loop vacuum integrals, namely the p parameter 
at (9(aa^) p, |^. The remaining cases have been considered in ^ ^, |1^, |1T|. Thus it is 
— at least in principle — possible to treat all problems where exactly one heavy mass is 
involved. 

In this paper we want to present the program package MATAD which was designed for 
the computation of MAssive TADpoles at one-, two- and three-loop order as pictured in 
Fig. |l|. Thereby each line may be massless or carry the mass M. In mathematical form 
the integrals to be solved by MATAD read 

d^p 1 



(27r)^ (|j2 + M2)' 
d^p d^k 



(27r)^ {27i)D (pI + MlY^{j)l + MlY-{j>l + M|)"3 




Figure 1: Prototype topologies for one-, two- and three-loop vacuum diagrams. The 
momentum pi flows through the line i as indicated by the arrow. Each line may either be 
massless or carry mass M and may be raised to an arbitrary integer power. 
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with Mj = or Mi = M. These expressions correspond to the diagrams in Fig. |l] where 
the momentum pi flows through the line i as indicated by the arrow, pi can be expressed 
as a linear combination of the loop momenta. However, these relations are in our case 
not of interest. Note that the integrals in Eq. (|l|) are defined in Euclidean space. 

The key idea for the computation of tadpole integrals is based on the integration-by- 
parts method [^ (see also Appendix ^. It can be used for the derivation of recurrence 
relations which relate vacuum integrals with different denominator structures. The proper 
use of the recurrence relations allows the reduction of an arbitrary integral to simple 
ones, which can be solved by successively using one- and two-loop formulae, and a linear 
combination of a few so-called master integrals. Only for them a hard calculation is 
necessary. In the case of three-loop tadpole diagrams nine master integrals are needed. 

At first sight the applications for vacuum integrals seem to be quite restricted. How- 
ever, for diagrams involving several mass scales, which follow a certain hierarchy, it is very 
often advantageous to apply an asymptotic expansion [jl2| which allows for a systematic 
expansion in the inverse heavy scale. Then the multi-scale integrals are expressed as prod- 
ucts of single scale ones. In ||13[| and |]n[ the rules for the so-called large-momentum and 
hard- mass procedure have been automated and computer programs, LMP [|13] and EXP |14|, 



have been developed. They generate for a given diagram all relevant sub-graphs together 
with the administrative files which govern the very calculation. 



LMP 



|I^ is written in PERL 

and can be applied to problems where one large momentum is involved. EXP ||I4|, writ- 
ten in C++, allows for a successive use of the large-momentum and hard-mass procedure 



and thus can deal with problems involving many scales. Both programs produce output 
which can be read into MATAD and MINCER. Thus the combination of both massive vacuum 
integrals and massless propagator-type diagrams is very powerful to attack problems in- 
volving several different mass scales. We want to mention that MATAD can be easily linked 
to a generator for Feynman diagrams. More details — in particular on the automation of 
the computation of Feynman diagrams — can be found in |jl5| . 



The outline of the paper is as follows: In Section |^ the structure of MATAD and the 
way it works is described. With the help of this section the reader should be able to 
use MATAD for his own problems. Deeper insight into some selected parts is provided in 
Section |. In Section H explicit examples are discussed and hints for the convenient us- 
age of MATAD are given. In Appendix ^ the ideas of the integration-by-parts method are 
reviewed. Appendix ^ lists all massive/massless combinations which are implemented 
into the topology files and in Appendix |C| the notation of the input and output is de- 
scribed. Furthermore the results for the master integrals are listed and the switches for 
the input-file are described. Appendix contains a list of all files of MATAD and, finally, 
in Appendix |^ the complete output of one of the considered examples is listed. 

2 Structure and mode of operation 

As MATAD is completely written in FORM [Q its installation reduces to copying the individual 
files into the corresponding directories. In the main directory the following files appear: 

form. set inc/ matadform pre/ problems/ 

The directories inc and pre contain the include-files and procedures, respectively. They 
are described in more detail in Section ^ and Appendix 0. matadform is a shell script 
which calls FORM in such a way that files from sub-directories can be included. It has to 
be adjusted by the user by simply specifying the corresponding paths. The file form. set 
contains FORM-specific settings which have to be adjusted according to the underlying 
platform. For details concerning the different switches we refer to the FORM manual 0]. 
The user-specific files are all contained in the folder problems. 

There are at least two files which should be provided by the user: main<prb> and 
<prb>.dia where <prb> stands for the name of the considered problem. The first one 
contains apart from some parameters essentially the information which diagram should be 
treated. Some explicit examples are given below. All the information about the diagrams, 
the projectors to be applied, etc. is contained in the file <prb> . dia. It is built up by FORM 
folds and splits into two parts so that its generic structure looks as follows 

*— #[ TREATO: 

[...] 

*— #] TREATO: 

*— #[ TREATl: 



*— #] TREATl: 

*— #[ TREAT2: 

[...] 

*— #] TREAT2: 

*— #[ TREATMAIN: 

[...] 

*— #] TREATMAIN: 

* 

* in the following list each diagrami is contained in a separate FORM fold 

* 

*— #[ dill: 
[ . . . diagrcun 1 . . . ] 
#define TOPOLOGY "XY" 
*— #] dill: 

*— #[ dll2: 
[ . . . diagrsun 2 . . . ] 
#define TOPOLOGY "XY" 
*— #] dll2: 

[...] 

The first part consists of tlie first four folds — tlie so-called special treat files. They 
provide the possibility to interact at different stages and thus influence the computation. 
Whereas TREATO, TREATl and TREAT2 are read before the recurrence relations are applied 
the content of TREATMAIN is read right before the results are stored to disk. The second 
part of <prb>.dia contains a list of all diagrams to be considered where each diagram is 
written in a separate fold. The name of these folds is arbitrary. 

Once these two files are set up the calculation is simply initiated by calling the program 
main<prb> and the following steps are performed. They are also illustrated in Fig. 0. 

1. Read global settings. They are partly contained in inc/main . gen and should not be 
modified. Others can be set by the user in the file main<prb>. They are described 



in Appendix C.4 



2. Read the input data for the diagram specified in main<prb> with the help 
of the variables PRB, FOLDER and DIAGRAM. The generic FORM command reads 
#include problems/ ' PRB V FOLDER ' .dia # 'DIAGRAM'. 

As a next step the file treat .pre is called and the following operations are performed. 

3. Insert Feynman rules for functions appearing in the input. In a first step the fermions 
(encoded in the functions S, SS, ... , cf . Appendix |(J.1D are resolved. It is important 
to do this before any contraction of indices is done. Then the propagators and 
vertices are treated. 
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read diagram 
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call treat .pre: 

• treat fermion lines 

• insert Feynnian rules 

• expand in small momenta 

• evaluate traces 

• perform Wick rotation 

• apply derivatives 

• read special treat files 



include input top file: 

• map expression to basic topologies 



include basic top file: 

• treat scalar products in numerator 

• apply recurrence relations 

• insert master integral 



TREATMAIW 



expand in e and store result to disk 



Figure 2: Flowchart illustrating the structure of MATAD. 



In the current version the QCD Feynnian rules are implemented (except the four- 
gluon vertex; see Appendix p.l|) . It is, however, straightforward to implement new 
vertices in the user-specific treat files. 

4. Apply projector. This should be done in one of the special treat files. The optimal 
position depends on the integrals to be computed. From now on only scalar integrals 
without any free indices are present. 



5. Expansion of the scalar denominators in the small quantities (mass and/or momen- 



turn). 

6. Perform traces. 

7. Do Wick rotation. This is done by multiplying each momentum by the imaginary 
unit (see also Appendix [C]). From now on the expression is defined in Euclidean 
space. 

8. Apply derivatives in order to factorize the external momentum. In this context see 
also the variables DALA12 and DALAQN in Appendix |C.4| . 

As there is the possibility to interact at three different places — after the fermions are 
treated (TREATO) and before and after the traces are performed (TREATl, respectively, 
TREAT2) the order of the commands may slightly be varied by the user. 

At this stage the scalar products in the numerator of the integrals should be formed 
by either only loop momenta or only external momenta (which then constitutes a trivial 
prefactor). In the denominator the (scalar) propagators may be raised to arbitrary power. 

The next steps constitute the main part of MATAD. 

9. Express the scalar products of the numerator in terms of the denominators. This 
produces a "1" in the numerator of the integrals. It might be that this step is very 
time and memory consuming. 

10. Apply recurrence relations to reduce the number of different integrals to simpler 
ones and to a small set of master integrals. 

11. Expand the result in e and store it in the directory 

problems/ 'PRB '/results/ 'DIAGRAM' .res under the name 'DIAGRAM'. 

An expansion in e is also done at various intermediate steps. Although poles of at most 
third order may appear for a three-loop vacuum integral terms up to order e^ have to 
be kept in the expansion as artificial poles may appear during the application of the 
recurrence relations (cf. step 10). 

Steps 9 and 10 constitute the central part of MATAD. They heavily depend on the loop- 
order and the topology which has to be specified apart from the very diagram in the 
folds dill, dll2, . . . (see above). Thus let us elaborate on this point in the following. In 
principle it suffices to define one input topology at one-, two- and three-loop order, where 
the number of internal lines amount to one, three and six, respectively (see Fig. |1|). If 
one allows each line to be massless and carry the mass M at the same time these three 
topologies are sufficient to cover all possible cases that can occur in the calculation of 
one-, two- and three-loop vacuum integrals. Note that a partial fractioning for terms like 



{V^f ij)^ + M2) 



b ' 



(2) 



where a and h are positive integers, leads to the same topologies with the only difference 
that now each line is either massless of massive. It is, however, not at all practical to 




Figure 3: The input of the diagram to be computed can be mapped to one of these 
topologies. The implemented massive/massless combinations can be found in Appendix p. 
The momentum p^ flows through line i as indicated by the arrow. 



rewrite the input to the notation of Fig. |1| before generating the file dia.<prb>. On 
the contrary it is advantageous to enlarge the input topologies. Currently the topologies 
shown in Fig. ^ are implemented in MATAD. The momentum pi flowing through line i 
can be expressed as a linear combination of the loop momenta. For our purpose these 
relations are, however, not of interest. The choice of the topologies was guided by the 
package MINCER [Q and for convenience the same notation concerning the definition of 
the momenta pi has been adopted. The implemented massive/massless combinations of 
each topology are listed in Appendix ^. After the declaration of the diagram in the folds 
dill, dll2, . . . the corresponding topology is specified via 

#define TOPOLOGY "XY" 

where XY corresponds to one of the topologies of Fig. ^ 

The notation concerning the momenta as introduced in Fig. |^ is quite convenient to 
be used for the input. However, the very recursion procedure is formulated for the three- 
loop topology of Fig. |l] where the lines are either massive or massless. This leads to 14 
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Figure 4: Basic three-loop topologies implemented into MATAD. They are reduced to 
either simple integrals or master integrals with the help of recurrence relations. 



different cases which are classified in P]. Thus before step 9 is performed the momenta are 
transformed from the notation of Fig. ^, which is used in the input, to the so-called basic 
topologies shown in Fig. ^. For them — after decomposing the scalar products in the 
numerator into parts of the denominator — the very recursion procedure is performed. 
Note that at this stage all propagators may be raised to an arbitrary integer power. 

Actually some of the three-loop diagrams (e.g. BN3) can be computed by the suc- 
cessive use of one- and two-loop procedures for massless propagator type diagrams or 
vacuum integrals, respectively. In such cases some of the (one- and two-loop) routines 
from MINCER [Q are used for parts of the computation. The corresponding procedures 
are listed in Appendix 0. For other cases (e.g. E4 or BN2) simple relations reduce one 
of the lines to zero and the resulting diagram can again be solved easily. Only for the 
cases D5, D4, DN, DM, E3, BN and BNl the recursion procedure has to be applied until 
one arrives at master integrals. They coincide with the corresponding diagrams of Fig. ^ 
where all denominators are raised to power one. Only the one pictured in Fig. ^, which 
results from BNl, is needed in addition. From this diagram even the 0{e) part is required. 
The analytic expressions are given in Appendix |C.3| . It should be mentioned that the 
recurrence relations for BM are quite involved. However, for this topology no difficult mas- 




Figure 5: Master diagram resulting from topology BNl. Here all propagators are raised 
to power one. Its constant part contains the expression S2 and its 0{e) part 0epS2 as 
listed in Appendix |C.3| . 






Figure 6: Simple integrals which occur as a result of the recurrence procedure next to 
the master integrals. All propagrators may be raised to an arbitrary power. 



ter integral is needed. Note that the basic topologies BMl and BM2 actually coincide with 
BN2 and BN3, respectively. For convenience partly different codes had been written which 
actually turned out to be quite useful while debugging MATAD. 

The recurrence relations leading to simple diagrams and master integrals for the three- 
loop topologies are derived in Refs. f^ and ^. Except the case where all six lines are 
massive all of them are implemented into MATAD. The reason why this case is missing is 
simply that it was not yet needed for practical calculations. All master integrals are listed 
in Appendix |C3| . 

In Fig. P the simple integrals are listed which also result from the recursion procedure. 
They can all be expressed in terms of F functions by the successive use of results for 
massless one-loop two-point functions, Pab{Q), in combination with one- (Va) and two- 
loop [Vabc) vacuum integrals. For convenience we want to list the explicit results for 
Pabil), Va and Vabc in Euclidean space: 
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^^2^D/2^a-b r(^ ^ ^ _ n/2)r{D/2 - a)r{D/2 - b) 

(Airf/' T{a)T{b)T{D - a - b) 



^" ' {2'Kf{p' + M^y (4vr)^/' r(a) 
r d^p d^k 1 



(27r)^ (27r)^ (p2 + m2)'^(F + M^fdp + k)y 
^^2^D-a-b-c p(^ _^ ^ ^ ^ _ ^)p(^ ^ ^ _ /)/2)r(6 + c - D/2)r(D/2 - c) 
(47r)'' r(a)r(6)r(a + 6 + 2c - D)T{D/2) " 

(3) 




Figure 7: Three-loop diagram of type LA. 



The vacuum integrals occurring at one- and two-loop level are quite simple. Actually 
most of them can be expressed in terms of F functions for arbitrary exponents of the 
propagators and no recursion relations are needed (cf. Eq. (^). Only for the two- loop 
integral in Fig. |l| where all three lines carry the mass M it is useful to implement simple 
recurrence relations which reduce the integrals to one master integral where all exponents 
are raised to the first power only. For three-loop calculations the result of this integral is 
needed up to 0{e) (Tlep, see Appendix |C3| ). 

For definiteness we want to consider an explicit example. Let us consider the three- 
loop diagram of Fig. ^ which we would like to expand up to fourth order in the external 
momentum, qi. For simplicity we neglect the tensor structure and consider only the scalar 
integral obtained in the case when the full line represents a massive particle and the curly 
line a massless one. In the directory problems/scalar/ one can find the following files 

mainscalar results/ scalar. dia 

where results/ is a directory to store the result of the diagram. The file scalar. dia 
looks as follows 
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*- 


-#[ TREATO: 


*- 


-#] TREATO: 


*- 


-#[ TREATl: 


*- 


-#] TREATl: 


*- 


-#[ TREAT2: 


*- 


-#] TREAT2: 


*- 


-#[ TREATMAIN: 


*- 


-#] TREATMAIN: 


*- 


-#[ scalar: 




M-4 




*slm 




*s2m 




*s3m 




*Dli(p4,ql) 




*Dh(p5,ql) 




*Dli(p6,ql) 




/p7.p7 




/p8.p8 



#defiiie TOPOLOGY "LA" 
* — #] scalar: 

In this case no special treat file is needed. The very diagram can be found in the fold 
scalar where the multiplication with M^ is done in order to end up with a dimensionless 
expression. The external momentum is routed through the lines 4, 5 and 6 as can be 
seen in the arguments of the function Dh. The other massive propagators are denoted by 
slm, s2m and s3m. The massless lines translate into the factors l/p7.p7 and l/p8.p8. 
For more details on the notation we refer to Appendix |C.1| . The file mainscalar looks as 
follows: 

#define PRB "scalar" 
#define DALAQN "ql" 
#define GAUGE "0" 
#define POWER "4" 
#define CUT "0" 
#define FOLDER "scalar" 
#define DIAGRAM "scalar" 
#- 
#include main. gen 

With the command 

#define POWER "4" 

we require an expansion up to fourth order in qi and 
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#define DALAQN "ql" 

effectively factors out the terms (qf)"'- The calculation is initiated with the command 

> matadf orm problems/scalar/mainscalar 

and after a few seconds the final result is displayed on the screen 

scalar = 

- 2 + 2*z3 - qi.Ql*z3*M'-2 + 227/216*Ql . Ql*M~-2 + l/2*qi . qi'2*z3*M'-4 
- 1876/3375*Ql.Ql-2*M--4; 

As expected it is finite and contains three terms in the expansion in qf/M'^. The result 
is stored to the directory problems/scalar/results/ and can be read with the FORM 
command load. Details on the notation of the output and the used conventions can be 
found in Appendix |C[ 

3 Some details on the internal structure 

Once the recurrence relations are implemented it is in principle possible to compute 
diagrams of arbitrary complexity. However, in practice one arrives quite soon at the 
limits set by the soft- or hardware. The algebra language FORM is designed to deal with 
large expressions. Still it happens quite easily that in internal steps the expressions exceed 
several Megabytes and even approach the order of a few Gigabyte. This significantly slows 
down the performance and it is advantageous to implement several tricks. Some of them 
are described in this section. 

• During the recurrence procedure it happens very often that whole blocks of com- 
mands have to be executed until the recursion has reached an end. Within FORM there 
are the commands 

repeat ; 
endrepeat ; 

which in principle allows for such a construction. However, in practice it is not possi- 
ble to use between repeat and endrepeat a command which forces FORM to combine 
identical terms like, e.g., .sort. For this reason a preprocessor variable, NOR (see also 
Appendix |(J.4|), has been introduced which in combination with the construction 



#do i=l, 'NOR' 
#enddo 

allows the use of . sort in intermediate steps of the recursion commands. 

• At this point we should also mention that the procedure ACCU has been adopted 
from MINCER Q]. It collects in the argument of the function acc() the polynomials in e 
and thus significantly reduces the number of terms. E.g., the expression 
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xl*x2 + 4*ep*xl*x2 + 12*ep~3*xl*x2 
transforms after #call ACCU{test} to 
accd + 4*ep + 12*ep~3)*xl*x2 

and instead of three only one term has to be treated in the following commands. 

• Very often it happens that propagators of the type 1/(M^ — {p + lYY have to be 
expanded in the momentum q and afterwards derivatives w.r.t. q are applied in order to 
factor out powers of q^. It turns out that it is very useful to expand in a first step only the 
part 2pq and keep the factors q^ unexpanded in the form 1/(M^ — (p^ + g^))". Thus less 
terms have to be considered while the derivatives are applied. Afterwards the expansion 
in q^ is performed. A related discussion can also be found in ||16|| . 



• The expansion of the scalar denominators in a small momentum is also very time 
consuming — especially for high values of 'POWER'. In order to do this in an effective way 
the variable poco, which is an abbreviation for "power counting", is defined via 

S poco(: 'POWER'); 

after the declaration of the preprocessor variable POWER. This definition ensures that the 
terms involving poco" with n > 'POWER' are automatically set to zero. Depending on the 
problem poco should also be considered in the special treat files. 

• The application of the recurrence relations can lead to spurious 1/e poles (cf. Ap- 
pendix ^ which are in general quite dangerous if an expansion in e is performed in 
intermediate steps. In MATAD at most three 1/e poles arise from the recurrence relations. 
Together with a possible 1/e^ term from the three-loop integrals an expansion up to order 
e^ has to be done in intermediate steps in order to get the correct constant term. 

• For quite a lot of applications the topology BN (cf. Fig ||) plays a crucial role. The 
original recurrence procedure for this toplogy to master integrals was proposed in 0: In 
a first step the exponent of three out of the four massive denominators are reduced to 
one. Then the exponents of the massless lines are reduced to zero. Finally the remaining 
line is treated and one arrives (apart from simple integrals) at an integral consisting of 
four massive lines connecting two vertices. It is connected to the corresponding master 
integral of Fig. ^ through a simple relation. 

The equation involved in the last recursion step is quite involved. It generates from 
each term more than ten terms at each call. Note that the exponent of the last mas- 
sive denominator gets increased by the proceeding steps. This enormously slows down 
the calculation — in some cases it makes it even impossible. The idea to circumvent 
this problem is based on the observation that the massless exponents can be reduced 
to zero even if none of the massive ones is reduced to one. The corresponding recur- 
rence relations are short and thus one arrives with only little effort at three-loop integrals 
with four massive lines, i?7v(0, 0, 71,3, 77,4, ^5, rig) (cf. Fig. ^ where all lines are massive and 
the exponents of the denominators are given by ^3, ^4, n^ and ng). These diagrams are 
now treated in the following way: Temporarily an external momentum is introduced which 
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flows through one of the lines. We choose line 3 for definiteness. In a second step the oper- 
ator Dg = djdq'^d/dq^ is applied and q is set to zero afterwards. This leads to an equation 
connecting -BAr(0, 0, ris, 77,4, 77,5, rig), -BAf(0, 0, ris — 1, ni.n^.no) and Bn{Q, 0, ^3 — 2, n^^.n^^no) 



1 
AM\n;~2){n^-l) 



BN{0,0,n3,n4,n5,n(i) = -ytttT. ^vTTI 7T°g^^(0' 0' ^3 - 2,^4,^5, ng) 



-2L) + 4(^3 - 1) , X / X 

H -^^^ — _ Ejv(0,0,n3- 1,^4,^5, ne), (4) 

which can be applied until 713 = 3. As the other indices are not affected the same procedure 
can be applied to them as well and one ends up with the integrals -BAr(0, 0, 1, 1, 1, 1), 
Bn{0, 0, 2, 1, 1, 1), Bn{0, 0, 2, 2, 1, 1), 5^(0, 0, 2, 2, 2, 1) and S^(0, 0, 2, 2, 2, 2). Both their 
values for g = and the result for the application of (Clg)" has to be known where the 
index n depends on how often Eq. (^ had to be applied. The overall number of the 
integrals needed is still small and for most practical applications well below 100. They 
can be computed once and for all using the method of Ref. and can be collected in a 
table. Currently the table contains all results up to n = 10 and partly for n = 11. In case 
the table is too small the original recurrence procedure |^ has to be used. This is done 
by defining the preprocessor variable BNRECOLD in the main file. 

4 Examples 

In this section we want to discuss some typical examples which can be treated with 
MATAD. In particular the content of the TREAT folds and the switches in main . <prb> shall 



be discussed in detail. In Section ^]T| two-loop corrections to the photon polarization 



function are considered and in Section [4.2| the calculation of three-loop vertex corrections 



contributing to the Higgs decay into gluons is discussed. As a last example we consider 
the computation of the fermion propagator in the limit of a small external momentum. 

4.1 Photon polarization function 

In this section only a calculation of two-loop diagrams is presented. However, we want 
to take this opportunity to show a concept which can also be used for the treatment of 
more complex problems. In particular it is shown how the complete calculation can be 
automated and one can ensure that all results are indeed up-to-date. 

To be precise, we want to consider the two-loop diagrams induced by a massive quark 
to the photon polarization function. The problem shall be called Pi. Then the file Pi .dia 
looks as follows: 

* 

* problem: Pi 
*— #[ TREATO: 
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#message project out transversal part 

multiply, (d_(mul,mu2)-ql(mul)*ql(mu2)/ql.ql)*deno(3,-2) : 
. sort 
*— #] TREATO: 

*— #[ TREATl: 
*— #] TREATl: 

*— #[ TREAT2: 
*— #] TREAT2: 

*— #[ TREATMAIN: 
*— #] TREATMAIN: 

* 

* 2-loop diagrams for problem Pi 



*— #[ d211: 

((-1) 

*Dg(nul,nu2,p5) 

*S(mul ,ql ,plm,nul ,ql ,p2m,mu2,p3m,nu2,p4m) 

*1); 

#define TOPOLOGY "Tl" 
*— #] d211: 

*— #[ d212: 

((-1) 

*Dg(nul,nu2,p4) 

*S(mul ,ql ,-p2m,mu2,plm,nu2,p3m,nul ,plm) 

*1); 

#define TOPOLOGY "T2" 
*— #] d212: 

*— #[ d213: 

((-1) 

*Dg(nul,nu2,p4) 

*S(nul ,p3m,nu2,plm,mu2,-ql ,-p2m,mul ,plm) 

*1); 

#define TOPOLOGY "T2" 
*— #] d213: 

The function deno(x,y) means l/(a; + ye) and can be used for denominators. TREATO 
contains the projector to the transversal part and the three contributing diagrams are 
named d211, d212 and d213. The external momentum is denoted by ql. The mass of the 
quarks, which in the final result appears as M, is introduced through adding the symbol 
m to the corresponding momentum which is defined through the topologies Tl and T2 in 
Fig. ^. We want to assume that qf <^ M^ and thus perform an expansion in gi. This is 
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achieved with the notation S(. . . ,ql,plm, . . .). For more details concerning the notion 
we refer to Appendix |C[ 
The file mainPi reads: 

#define PRB "Pi" 
#define PRDBLEMO "1" 
#define DALAQN "ql" 
#define GAUGE "xi" 
#define POWER "4" 
#define CUT "1" 
#define FOLDER "Pi" 
***#define DIAGRAM "d211" 
#- 
#include main. gen 

For most of the specified variables we refer to Appendix |C4| . We only want to mention that 
a general gauge parameter xi is chosen with the command #def ine GAUGE "xi". This 
allows for an explicit check that the final result, i.e. the sum of the three diagrams is gauge 
parameter independent. The definition #def ine POWER "4" requests for an expansion up 
to fourth order in ql. The definition of the variable DIAGRAM is commented as it will be 
defined during the call of mainPi (see below). 

The computation of each diagram could be started separately and the results could be 
summed at the end. Instead we want to take the opportunity to present a method which 
is unavoidable for problems where a large number of diagrams contribute. In the following 
we want to present two more files which can easily be adopted to other problems. The 
first one, makePi, is a so-called GNU make file and could look as follows: 

SHELL = /bin/sh 

DIA2 = \ 

problems/Pi/results/d211 . res\ 
problems/Pi/results/d212 . res\ 
problems/Pi/results/d213 . res 

problems/Pi/results/Pi . 2 . res : $ (DIA2) 

matadf orm problems/Pi/comPi > problems/Pi/log/comPi . log 

ii = $(notdir $(basename $@)) 

$(DIA2): problems/Pi/mainPi 

if [ -f problems/Pi/results/$(ii) .res ]; \ 

then rm problems/Pi/results/$(ii) .res; fi 
time matadf orm -d DIAGRAM=$(ii) \ 

problems/Pi/mainPi > problems/Pi/log/$(ii) . log; 
if [ -f problems/Pi/results/$(ii) .res ]; \ 
then rm problems/Pi/log/$(ii) . log; fi 



For details concerning the individual commands we refer to the literature [17|. For us 



it is only important that at the beginning the diagrams we want to compute are listed. 
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Furthermore, after the hne 'problems/Pi/results/Pi . 2 . res : $ (DIA2) ' the command is 
given which specifies what shall be done once the computation of the individual diagrams 
is finished: The diagrams are summed with the help of the program comPi 



* 

* comPi 

* 

#- 

#include declare .matad 

#define PRB "Pi" 

.global 

#do i=l,3 

load problems/ ' PRB ' /results/d21 ' i ' . res ; 
#enddo 

g res 'PRE '2 = 
#do 1=1,3 
+ d21'i' 
#enddo 



b ep; 
print ; 
. store 

#+ 

save problems/ ' PRB ' /results/res ' PRB ' .2. res res ' PRB ' 2 ; 

. end 

To initiate the calculation one simply has to specify the command 

> make -f problems/Pi/makePi 



where mcike has to call the GNU version fl^ of the make command. The final result is 



stored in the file problems/Pi/results/resPi . 2 . res. It reads 

resPi2 = 

+ ep'-l * ( - 6*Q1.Q1 + 8/5*Ql.Ql~2*M~-2 ) 

+ ep * ( - 35/6*Ql.Ql - 6*Ql.Ql*z2 + 8/5*Ql . Ql~2*M~-2*z2 + 3116/1215* 
qi.qi'2*M'-2 ) 

+ 13/3*Q1.Q1 - 128/405*Ql.Ql-2*M--2; 

Indeed, as expected, there is no trace of the gauge parameter xi. 
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4.2 Higgs decay into two gluons 

Party for convenience and partly for historical reasons the notation of the input topologies 
in Fig. ^ is closely connected to two-point functions. However, MATAD only deals with 
vacuum diagrams independent of the number of external legs. In this section we want to 
show that also problems involving at first sight three-point functions can be approached 
using MATAD. 

Let us consider QCD corrections to the decay of the Standard Model Higgs boson 
into two gluons. It is convenient to construct an effective theory where the top quark is 
integrated out. Details on the theoretical background can be found in [|l^. Here it shall 
only be mentioned that the coefficient function which contains the dependence on the 
mass of the top quark, Mj, can be computed from triangle diagrams as pictured in Fig. ^. 
According to their Lorentz structure the result can be written as follows 

K{M,) {qU^ - q.q^gn , (5) 

where gi and g2 are the momenta of the gluons with polarization vectors e'^(q'i) and e'^(q'2)- 
Thus the vertex diagrams have to be expanded up to linear order both in gi and q2 and 
an appropriate projector has to be applied in order to get K{Mt). 

J)J)J)59J)Jb 

Figure 8: Sample diagram contributing to the decay of the Higgs boson. Solid and looped 
lines represent quarks and gluons, respectively. 



The file containing the diagrams could look as follows 
* 

* problem: hgg 
* 

*— #[ TREATO: 
multiply, ( 

a*deno (2 , -2) * (ql . q2*d_ (mu ,nu) -q2 (nu) *ql (mu) -q2 (mu) *ql (nu) ) 
+b*deno (2,-2)* (-ql . q2*d_ (mu , nu) + (3-2*ep) *q2 (nu) *ql (mu) +q2 (mu) *ql (nu) ) 

); 

. sort 

*— #] TREATO: 

*— #[ TREATl: 
*— #] TREATl: 
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*— #[ TREAT2: 
*— #] TREAT2: 

*— #[ TREATMAIN: 
*— #] TREATMAIN: 

*— #[ d31335: 
((-1) 
*M 

*Dg(iiul,nu2,pl) 
*Dg(iiu7,nu8,-p4) 
*Dg (nuS , iiu4 , ql , -pi ) 
*Dg (iiu5 , iiu6 , q2 , p 1 ) 

* S ( -ql , -p3m , nu? , -ql , p5m , nu4 , -p2m , nu6 , q2 , p5m , nu8 , q2 , -p3m) 
*V3g(mu,ql ,iiul ,-pl ,nu3,pl-ql) 
*V3g (nu , q2 , iiu2 , pi , iiu5 , -pl-q2) 
*1); 

#define TOPOLOGY "04" 
*— #] d31335: 

where the diagram d31335 corresponds to the one shown in Fig. |^. The fold TREATO 
contains (up to an overall factor {qi-q2)~'^) the projector on the coefficients in front of the 
structures g^'^ and q^ql^ of Eq. (|^). They are marked by the symbols a and b, respectively. 
Thus the transversality of Eq. (|^) can be explicitly checked in the sum of all contributing 
diagrams (the result of a single diagram does in general not have a transverse structure). 
The corresponding main-file is very similar to the one listed in Section [4.1| . The only 
difference (apart from replacing Pi by hgg) is that the commands 

#define DALAQN "ql" 
#define GAUGE "xi" 
#define POWER "4" 
#define CUT "1" 

should be replaced by 

#define DALA12 "1" 
#define GAUGE "0" 
#define POWER "2" 
#define CUT "0" 

The third line ensures that an expansion of the integrand up to second order in the 
external momenta is performed and the first one sets qf and gg to zero and factors out 
the scalar product qiq2- #def ine CUT "0" sets e to zero in the final result as the terms 
of 0{e) are anyway not computed completely at three- loop order. In this example we 
choose Feynman gauge which is achieved with #def ine GAUGE "0". 

After calling MAT AD it takes of the order of a minute to obtain the result: 
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d31335 = 

+ ep--2 * ( 40*qi.q2*M-2*a + 344/9*Ql . Q2-2*a - 232/9*qi . q2-2*b ) 

+ ep--l * ( - 308/3*Ql.Q2*M-2*a - 3530/27*Ql . Q2-2*a + 1786/27*qi . q2-2* 
b ) 

+ 60*qi.q2*M'2*z2*a + 734/3*qi . q2*M'2*a - 1936/9*qi . q2'2*z3*a + 1136/9* 
qi.q2-2*z3*b + 172/3*qi.q2-2*z2*a - 116/3*Q1 . Q2-2*z2*b + 46817/81* 
qi.q2-2*a - 26239/81*Ql.Q2-2*b; 

Note that the terms proportional to Q1.Q2*M~2 cancel after adding all contributing dia- 
grams. 

4.3 Fermion propagator 

In this example we compute the small-momentum expansion of a three- loop diagram 
which contributes to the fermion propagator. Recently these kind of diagrams have been 
considered in different kinematical regions in order to obtain the three-loop relation be- 



tween the MS and on-shell quark mass [1^, [Tl|]. This subsection contains all relevant input 



information whereas the complete output is given in Appendix |^. 

The fermion self energy, S(g), can be decomposed into a scalar and vector part 

S(g) = MS5(g') + ^Sv(g'), (6) 

where q is the external momentum and M is the mass of the quark. We are interested in 
the computation of the scalar functions S5 and Sy. 
The file mainf p looks as follows 

#define PRB "fp" 
#define PROBLEMO "1" 
#define DALAqN "ql" 
#define GAUGE "0" 
#define POWER "1" 
#define CUT "0" 
#define FOLDER "fp" 
#define DIAGRAM "d3179" 
#- 
#include main. gen 

The variable POWER is defined in such a way that an expansion up to third order is 
performed. This leads to the constant and order g^/M^ terms for the functions S5 and 
Sy as S(g) itself has mass dimension one. The corresponding projectors and the diagram 
to be computed look as follows (fp.dia): 

* 

* problem: fp 
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*— #[ TREATO: 

multiply, 1/4*(1/M + a + g_(l ,ql)/ql .ql) ; 
*— #] TREATO: 

*— #[ TREATl: 
*— #] TREATl: 

*— #[ TREAT2: 
*— #] TREAT2: 

*— #[ TREATMAIN: 
*— #] TREATMAIN: 

*— #[ d3179: 
((1) 

*Dg(iiu3,iiu4,p5) 
*Dg (iiu5 , iiu6 , -p6) 
*Dg (nul , iiu2 , -ql , -pi) 

*S (nu2 , -plm , nu5 , -p4m, nu4 , -pSm ,nu6 , -p2ni , nu3 , -plm , nul ) 
*1); 

Sdefine TOPOLOGY "04" 
*— #] d3179: 

Note that in the fold TREATO the vector part gets multiphed by a in order to distinguish 
Sy from Us in the final result. The input for the diagram corresponds to the one pictured 
in Fig. ^. The complete output appearing on the screen and the result can be found in 
the Appendix ^ 



^% 







Figure 9: Sample diagram contributing to the fermion propagator. Solid and looped lines 
represent quarks and gluons, respectively. 



Acknowledgments 

It is a pleasure for me to thank K.G. Chetyrkin, R. Harlander, J.H. Kiihn and T. Seiden- 
sticker for many discussions and useful advice within the recent years when MATAD was 
developed. This work was supported in part by DFG under Contract Ku 502/8-1 {DFG- 
Forschergruppe "Quantenfeldtheorie, Computeralgebra und Monte- Carlo- Simulationen") 

22 



and by SUN Microsystems through Academic Equipment Grant No. 14WU0148. I am 
grateful to the department of Theoretical Particle Physics of the University of Karlsruhe 
for the pleasant atmosphere during a visit when a major part of this project was carried 
out. 



Appendix 

A Integration-by-parts 

The method of integration-by-parts plays a fundamental role in the computation of multi- 
loop diagrams . For this reason we want to demonstrate its underlying idea by consid- 
ering a typical example. 

The integration-by-parts algorithm uses the fact that the D-dimensional integral over 
a total derivative is equal to zero: 



/d°p^/(p,...) = 



dp- 



(7) 



By explicitly performing the differentiations one obtains recurrence relations connecting 
Feynman integrals of different complexity. The proper combination of different recurrence 
relations allows any Feynman integral (at least single-scale ones) to be reduced to a small 
set of so-called master integrals. The latter ones have to be evaluated only once and for 
all, either analytically or numerically. 




Figure 10: Two-loop master diagram. The arrows denote the direction of momentum 
flow. 



Let us for definiteness consider the scalar two-loop diagram of Fig. ^. The corre- 
sponding Feynman integral shall be denoted by 



J(ni,...,n5) 



d> d^A; 



Dl 



1 



{2n)D (27r)^ {pi + mf)-! • • • (pi + mlY^ ' 
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where pi, . . . ,p5 are combinations of the loop momenta p, k and the external momentum 
q (we work in Euclidean space here), rii, . . . ,11^ are called the indices of the integral. 
Consider the sub-loop defined by the lines 2, 3 and 5, and take its loop momentum to 
he p = p^. If we then apply the operator (d/dp^) ■ p^ to the integrand of J, we obtain a 
relation of the form (0) , where 



(pi + mi)"5(pi + mi)"2(pi + 



m 



ns 



Performing the differentiation and using momentum conservation at each vertex one de- 
rives the following equation: 



-■n-sS (5 —4 + ?Ti4 — rrig — r/igj — ^22 (5 —1 + rrii — rri;^ — 1712 

+ D — 2n^ — n^ — n2 + 2n^m\h'^ /(rii, . . . , 77,5) = 0, (10) 

where the operators 1^,2^,... are used in order to raise and lower the indices: 
I^J(. . . , nj, . . .) = /(..., ?2.j ± 1, .. .). In Eq. (p!0D , generally referred to as the triangle 
rule, it is understood that the operators to the left of /(ni,...,n5) are applied before 
integration. If the condition m^ = 0, m^ = m^ and rui = m2 holds, increasing one index 
always means to reduce another one. Therefore this recurrence relation may be used to 
shift the indices ni, n^ or n^ to zero which leads to much simpler integrals. 

The triangle rule constitutes an important building block for the general recurrence 
relations. The strategy is to combine several independent equations of the kind (|10[) in 
order to arrive at relations connecting one complicated integral to a set of simpler ones. 
For example, while the direct evaluation of even the completely massless case for the 
diagram in Fig. |10| is non-trivial, application of the triangle rule (|10D leads to 



[ni,...,nr,} 



D - 277-5 - n2 - ns 



^22+ (5- - 1-) + 7138+ (5- - 4-) ] /(ni 



,,n5J. 

(11) 
Repeated application of this equation reduces one of the indices ni, 724 or n^ to zero. For 

example, for the simplest case {rii = n2 = . . . = rzs = 1) one obtains the equation pictured 

in Fig. |ll|: The non-trivial diagram on the l.h.s. is expressed as a sum of two quite simple 

integrals which can be solved by applying the one-loop formula. This example also shows 

a possible trap of the integration-by-parts technique. In general its application introduces 

artificial 1/e poles which cancel only after combining all terms. They require the expansion 

of the individual terms up to sufficiently high powers in e in order to obtain, for example, 

the finite part of the original diagram. This point must carefully be respected in computer 

realizations of the integration-by-parts algorithm: One must not cut the series at too low 

powers because then the result goes wrong; keeping too many terms, on the other hand, 

may intolerably slow down the performance. 

In our example, the l.h.s. in Fig. |l^ is finite, each term on the r.h.s., however, develops 

l/e^ poles. The first three orders in the expansion for e — >■ cancel, and the 0{e) term 

of the square bracket, together with the 1/e in front of it, leads to the well-known result 

(omitting factors I/IGtt^): /(1, 1, 1, 1, 1) = 6C(3)/g^, where q is the external momentum. 
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0-<x> 



Figure 11: Symbolic equation resulting from Eq. (|Tl|) applied to the diagram /(1, 1, 1, 1, 1). 
The dot indicates that the respective denominator appears twice. 



In general, the successive application of recurrence relations generates a huge number 
of terms out of a single diagram. Therefore, a calculation carried out by hand becomes 
very tedious and the use of computer algebra is essential. 

B The topology files 

This part of the appendix provides a complete list of those massive/massless combinations 
which are implemented in the topology files inc/TOPOLOGY/topXY. In the following tables 
for all three-loop topologies of Fig. ^ the lines are listed which have to be massive. All other 
lines may be massless or absent. In some cases some of the lines have to be completely 
absent, i.e. it is even not allowed for them to be massless. This is explicitly specified in 
the columns "absent". 



LI 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1 


— 


2 


— 


1,2 


— 






Tl 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2 
2,4 




1,2,4 
1 




1,2,3,4 
2 




1,4 
1,2,5 




T2 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,3 
3,4 




2 
3 




1,2,3 
1,2 




2,4 




BE 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,3,4,5 


— 


1,2,3 


— 


4,5 


— 






BU 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,3,6,7 


— 


4,5,6,7 


— 
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LA 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,3,4,5,6 
1,2,3 


I 


1,3,4,6,7,8 
1,3,7,8 


4,5,6 


1,6,7 
1,7 


4,5,6 


3,4,8 
3,8 


4,5,6 


NO 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,3,4,5,6 


— 


1,2,3 


— 










N2 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,3,4,5,6 


— 


1,2,3,4 


— 


5,6 


— 


2,3 


— 


N3 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,3,4,5 
3,4,5 


I 


2 

4,5 


I 


1,3,4,5 
1,3 


I 


3,4 

1,5 


I 


01 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,3,4,6,7 

1,2 

5 

1,5 

1,6,7 


— 


1,2,3,4 

1,5,6 

5,7 

1 

1,2,7 


3,4 
3,4 


6,7 
1,5,7 
1,5,7 
1,7 


3,4 


1,2,6,7 
5,6 
7 
6 


— 


02 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,3,4,7 


— 


3,4,7 


— 


1,2 


— 


7 


— 
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04 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,6,7 


— 


1,2,6 


— 


7 


— 


1,3,5,6 


7 


1,3,5 


— 


3,5 


— 


1,5 


— 


5 


— 


1,3,5,7 


— 


3,5,7 


— 


1,2 


— 


1 


— 


2 


— 


1,6 


— 


2,6 


— 


1,3,7 


6 


3,7 


6 


2,6 


— 


1,3,5,7 


6 


1,5,7 


6 


1,2,7 


6 


2,3 


7 


1,2,3 


7 


3,4 


7 


3,4,5 


7 


1,3,4 


7 


1,2,5 


^7 

1 


1,3,4,5,6 


7 


1,4,5,6 


7 


1,3,4,5 


7 


1,4,5 


^7 

1 


1,4 


7 


3,4,6 


7 


1,2,3,4 


7 


2,3,4 


^7 

1 


1,2,3,4,5 


7 


1,2,4,5 


7 


1,3,4,6 


7 


1,2,3,5 


7 


2,3,5 


7 


3,4,5,6 


7 


4,5,6 


7 


2,3,4,5 


7 


2,4,5 


7 


2,4 


7 


1,2,4 


7 


1,5 


7 


1,3 


7 


3 


7 


2,5 


7 


1,2 


7 


4,5 


7 


4 


7 


1,2,4,5,6 


7 


2,4,5,6 


7 


6 


7 


4,6 


7 


3,5,6 


7 


3,6 


7 


1,5,6 


7 


2,5,6 


7 


1,3,6 


7 


1,2,5,6 


7 


1,2,4,6 


7 


2,4,6 


7 


5,6 


7 


1,4,6 


7 






Yl 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,3,6 


— 


1,3 


— 


3 


— 


3,5,6 


— 


Y2 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,3,5 


— 


2 


— 


1,3,5 


— 


3 


— 


5 


— 


1,3 


— 


1,5 


— 






Y3 


massive 


absent 


massive 


absent 


massive 


absent 


massive 


absent 


1,2,3,5 


— 


1,3,5 


— 


2 


— 


5 


— 


3,5 


— 


1,3 


— 











C Some details on the notation 

C.l Notation of the input 

In this subsection we describe the notation to be used for the input. 

On one side of the following expressions the FORM notation is used whereas the other 
side displays the corresponding mathematical terms. It is always assumed that Pi {i = 
1, . . . , 9) is a loop momentum and Qi {i = 1, . . . , 3) is a small momentum in which an 
expansion is performed. 

If scalar integrals are to be computed with MATAD the following notation has to be 
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used: 

— > 1/pl.pl, l/p2.p2,... 

-^ Dl(pl,ql),... 

— > slm, s2m, ... 
> Dh(pl,ql),... 

Fermions are treated with the help of the (non-commutative) function S and the function 
Dg can be used for the gluon propagator: 



massless hues: 


1 1 


massive hues: 


-iPi+giy^-- 


M'^-p^ ' M2-p^ ' • 




M2-(pi+gi)2' • • • 



S(mu) 




= 7^ 


S(pl) 


SCplm) 




M2-p2 ) 


S(ql,pl) 


S(ql,plm) 
Dg(mu,nu, 


Pl) 


-«2 


Dg(mu,nu 



-(Pi+9i)^ ' 



'' ^ -(Pl + gl)2 

-(pi+gi)^ 



(13) 



Instead of S also SS, SSS or SSSS may be used which is useful in case more than one 
fermion line is involved. The vertices between three gluons and gluon and ghosts are 
implemented through the following functions (where the colour factors are not taken into 
account): 



V3g(il,pl,i2,p2,i3,p3) = {P2 - PiY' g'''' + {P3 - P2y' 9'''' + {Pi - PsT' g'"'' 
Vgh(il,pl) = -pi^ 



(14) 



If small momenta are present they can simply be added to pi, p2 or p^. In analogy to the 
expansion in a small momentum it is possible to introduce functions where an expansion 
in small masses can be performed. 

In the above expressions instead of ql also q2 or qS may be chosen as momenta in 
which an expansion is done. For the loop momenta pl, . . . , p9 is allowed. 

For the Feynman rules given above the user has to check the consistency with own sign 
conventions. In particular all unnecessary factors like, e.g., the strong coupling constant 
are omitted. Also the imaginary unit, i, is suppressed. Furthermore the Feynman rules 
are chosen in such a way that the vertices are proportional to i and all propagators are 
proportional to 1/i which leads to convenient cancellations. Thus each diagram can be 
written as 

i \ — jy X [expression formed by the Feynman rules of (13) and (|14D ] , (15) 

where / is the number of loops. The factor 1/i in front of the integration momenta 
disappears after the Wick rotation. The overall factor i, which occurs for all diagrams 
and which is independent of the number of loops, is omitted. 
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C.2 Notation and conventions of the output 

At this point some words about the notation of the output of MATAD are in order. As 
already mentioned, the final expressions are expanded in e where — following general 
MS conventions — the factors 7^; and In 47r have been dropped. Furthermore the factors 
(l/167r^)' and {1/M'^'^Y where / is the number of loops are omitted. Concerning the 
remaining symbols the translation is given in the following table: 



Euclidean external momenta 
C2, Cs) ■ ■ ■ (Rieman zeta function) 
mass appearing in the integrals 
parametrization of the finite parts 
of the master integrals 



ep 

Q, Ql, Q2, ... 

z2, z3, . . . 

M 



D5, D4, ... (see Appendix ^) 

Note that all momenta which occur in the output have to be interpreted in Euclidean 
space. 



C.3 Master integrals 

MATAD does not insert the complete finite parts of the master integrals automatically. 
Instead they are parameterized by the symbols D5, D4, DM, DN, B4, ... where the notation 
is adopted from the one introduced in Fig. ^ In the following list one can find the 
corresponding analytical expressions §, ||, |10|, |ll], ^ ^ 



D6 
D5 

D4 
D3 
DM 
DN 
B4 

E3 



6C3-17C4-4C2ln'2 + -ln^2 



^ 469 ^ 8 



Clo - 



TT 



16Li4 

16 E ■ 

m>n>0 



-4 



CI, 



TT 



l2 



-l)'"cos(27m/3) 



-8.2168598175087380629133983386010858249695, 



6C3 



77 
12 



C4 



6C3 - ^C4 



6C3 



11. 



CI2 



CI2 hr 



vr 



TT 



Ghi-r 



TT 



1 

21 



6C3 - 4C2 In' 2 + - In^ 2 - y C4 + 16Li4 (^- 

2 13 /I \ 

-4C2 In' 2 + - In^ 2 - -C4 + 16Li4 (^- j , 

139 7ry31n^3 177r3v^ 21 1^ 

C2 + -Cs 

3 8 72 2 ^ 3^ 



^Note that there is a misprint in Eq. (22) of Ref. |1^: all four appearing Clausen functions must be 
squared and the coefficient of ^4 in the second to last equation must read "—22" instead of "—31/2". 
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TT 



-iir/a ' 



V^ 



S2 



0epS2 



+ 10V3Cl2f-) -6v^Im Lis 

763 gvrVs 111^3 SStt^v^ 195, 15, 57, 
32 16 48 16 ^^ 4 ^^ 16^^ 



Tlep 



45 vrv^ 111^3 
~Y 8 



Lis 



-i7r/6 ^ 



v^ 



^4c. + C, + 6v^Cb(^)-6v^I,„ 



Lis 



-in/d^ 



V3 



(16) 



with 012(2;) = Im[Li2(e*^)]. Li2, Lis and Li4 are the Di-, Tri- and Quadrilogarithm, 
respectively. S2 both appears in the diagram of Fig. |^ and the two-loop diagram of Fig. ^ 
where all three lines have the same mass. Their O^e) parts, which can contribute to the 
finite part of the three-loop results, contain 0epS2 and Tlep, respectively. We should 
mention that those expressions of Eq. (|16D which coincide with an existing topology (e.g. 
D5 <-s> topD5) indeed agree with the finite part of the corresponding master integral. On 
the other hand B4, e.g., comprises the complicated parts of the finite part of the master 
integral corresponding to topBN, however, not the complete one. D6 is not yet implemented 
into MATAD and listed for completeness only. 

C.4 Parameters and switches in main.<prb> 

In the following a brief description of the switches in the file main<prb> is given. In all 
cases the FORM syntax reads #def ine <VAR> "<VAL>" where <VAR> is the variable to be 
defined and <VAL> the corresponding value. 

It is required to define at least the following five variables: 

DIAGRAM: name of the diagram to be computed. 

FOLDER: the file containing the special treat files and the diagrams is called 
'FOLDER' .dia. 

GAUGE: determines the choice for the gauge parameter. The variable ^ as defined 
in the function Dg in Eq. (|T3|) is replaced by the value of GAUGE. In particular 
#define GAUGE "0" corresponds to Feynman gauge and with #def ine GAUGE "xi" 
the calculation is performed for general gauge parameter. 

POWER : determines the depth of the expansion in the small quantities. 

PRB: name of problem. PRE corresponds to the name of the fold where the problem- 
dependent files can be found. 
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The definition of tlie remaining variables is optional: 

BNRECOLD : if this variable is defined the original recurrence procedure of Ref. [^ is used 
for the topology BN. 

CUT: determines the depth of the expansion in e of the final result. At one-, two- and 
three-loop order at most the terms of order e^, e^, respectively, e° are reliable. The 
default value is "2" . 



DALA12: expansion in qiq2- If this variable is set ql and g| are set to zero and pow- 
ers in qiq2 are factored out. Currently only the expansion up to order {q\q2Y is 
implemented. Note that only positive powers of qiq2 can be treated. 

DALAQN: apply d'Alembert operator, Dg = d/dq^d/dq^j_, in order to factor out powers in 
q^. In the argument of DALAQN the momentum is specified with respect to which the 
derivatives are performed. 

NOR: is an abbreviation for Number Of Recursions. As the naive use of the repeat- 
endrepeat construction significantly slows down the performance the most compli- 
cated procedures are buffered by a #do-#enddo construction. 'NOR' constitutes the 
upper bound of the do-loop. The default value is "10". 

PROBLEMO/l/2/MAIN: If one of these variables is set a special treat file is read at the 
corresponding position. The value has to agree with the one defined in the file 
<prb>.dia. 

TIME: If this variable is set the statistics is printed at various steps of the calculation. 

There are more switches in inc/main.gen. However, they should not be modified as 
most of them are in an experimental stage and not sufficiently tested. 

D List of files 

In this appendix we provide a list of all files belonging to MATAD. The directory inc con- 
tains essentially the include-files. In particular some of the procedures are collected in 
files which are included at the very beginning. The tables for the topology BN are located 
in inc/TABLEDAL and inc/TABLEREC. The directories inc/TREAT and inc/TOPOLOGY es- 
sentially contain the files treating the individual (input and basic) topologies and in the 
folder pre some auxiliary procedures are collected. 

form. set inc/ matadform pre/ problems/ 

inc : 

TABLEDAL/ bnm2m expandDr nomBM reduceBM_2 tblBN 

TABLEREC/ declare. matad expepgaim nomBN reduceBN 

TOPOLOGY/ denoexp expnomdeno nomdecomBN reduceBNl 
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TREAT/ expandBN 
bnbm2prc expandBNM 
bnbmprc expandBNM_2 

inc /TABLED AL: 
BNd.tbl BNdO.tbl 



mam. gen 
mat ad. info 
matminprc 



recursion_2 

redcut 

redcutnomdeno 



reduceBN_2 
symmetryBM 
symmetryBN 



inc/TABLEREC: 

BN.tbl BNnl.tbl BNnln2.tbl BNn2.tbl 



inc/TOPOLOGY: 

topBE topBN 

topBM topBNl 

topBMl topBN2 

topBM2 topBN3 

topBM_2 topBN_2 topDN topMl topN2 top04 topY2 



topBU 


topE3 


topM2 


topNS 


topD4.add topY3 


topD4 


topE4 


topM3 


topNO 


topTl topempty 


topD5 


topLl 


topM4 


topOl 


topT2 


topDM 


topLA 


topM5 


top02 


topYl 



inc /TREAT: 
treat. dalal2 
treat . dalaav 
treatbm 
treatbml 



treatbm2 treatbn2 treatdS treate4 treatm4 

treatbm_2 treatbnS treatdm treatml treatmS 

treatbn treatbn_2 treatdn treatm2 treatnl 

treatbnl treatd4 treateS treatmS treattl 



pre : 

aver. pre cutep.prc dalaqn.prc difvecsc.prc solveS.pre treat. pre 

aver 1. pre dalal2sc.prc dalasc.pre pochtabl .pre tabBN.prc 



Files taken over from MINCER 

one. pre simplify. pre finish. pre tabtwo.pre 
aeeu.pre dotwo.pre newtwo.pre two. pre 
triangl2.pre triangle . pre poehtabl.pre 



E Fermion propagator: output 



In this appendix we present the complete output of the example discussed in Section ^]3 
Calling MATAD 

> matadf orm problems/f p/mainfp 

leads to 



FORM version 2.3 Apr 24 1997 



* mainfp 

* 

#define PRB "fp" 

#define PROBLEMO "1" 
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#define DALAQN "ql" 
#define GAUGE "0" 
#define POWER "1" 
#define CUT "0" 
#define FOLDER "fp" 
#define DIAGRAM "d3179" 
#- 
*~~ MATAD — computation of MAssive TADpoles 
*~~ read generic main file 
*~~ read diagrsun 
G dia= 
#include problems/ 'PRB'/ 'FOLDER' .dia # 'DIAGRAM' 

((1) 

*Dg(nu3,nu4,p5) 

*Dg (nu5 , nu6 , -p6) 

*Dg (nul , nu2 , -ql , -pi ) 

*S (nu2 , -plm , nu5 , -p4m , nu4 , -p3m , nu6 , -p2m , nu3 , -plm , nul ) 

*1); 

#define TOPOLOGY "04" 

*— #] d3179: 

#- 

*~~ Treat the traces 

*~~ Include special treat-file 

*~~ Feynman rules for vertices and propagators: 
*~~ gluon-ghost-ghost-vertex 
*~~ 3-gluon-vertex 
*~~ gluon propagator 
*~~ ghost propagator 

*~~ expand denominators 

*~~ Dh 

*~~ Dl 

*~~ 1 

*~~ 1 

*~~ Change notation to pl,p2, . . . 

*~~ Trace 1 

*~~ Trace 2 

*~~ Trace 3 

*~~ Trace 4 

*~~ treat DL(x) 

*~~ Do Wick-rotation 

*~~ Apply d Alembertian w.r.t. ql 

*~~ average done 

*~~ Expand Dr(p,q) 

*~~ 1 

*~~ q_i -> Q_i 

*~~ include TOPOLOGY-f ile 

*~~ this is top04 

*~~ Recursion of type d5 

*~~ this is topD5 

*~~ numerator 

*~~ do recursion 
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*~~ - done 

*~~ Recursion of type d4 

*~~ this is topD4 

*~~ numerator 

*~~ do recursion 

*~~ - done 

*~~ Recursion of type dm 

*~~ this is topDM 

*~~ numerator 

*~~ do recursion 

*~~ - done 

*~~ Recursion of type dn 

*~~ this is topDN 

*~~ numerator 

*~~ do recursion 

*~~ - done 

*~~ Recursion of type e4 

*~~ this is topE4 

*~~ numerator 

*~~ do recursion 

*~~ - done 

*~~ Recursion of type e3 

*~~ this is topE3 

*~~ numerator 

*~~ do recursion 

*~~ - done 

*~~ Recursion of type bn_2 

*~~ this is topBN_2 

*~~ numerator 

*~~ do recursion 

*~~ Use table for BN 

*~~ - done 

*~~ Recursion of type bnl 

*~~ this is topBNl 

*~~ numerator 

*~~ do recursion 

*~~ - done 

*~~ Recursion of type bn2 

*~~ this is topBN2 

*~~ numerator 

*~~ do recursion 

*~~ - done 

*~~ Recursion of type bn3 

*~~ this is topBN2 

*~~ numerator 

*~~ do recursion 

*~~ - done 

*~~ Recursion of type bm_2 

*~~ this is topBM_2 

*~~ numerator 

*~~ do recursion 

*~~ - done 
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*~~ 


Recursion of type 


bml 




*~~ 


this is topBMl 






*~~ 


numerator 






*~~ 


do recursion 






*~~ 


- done 






*~~ 


Recursion of type 


bm2 




*~~ 


this is topBM2 






*~~ 


numerator 






*~~ 


do recursion 






*~~ 


- done 






*~~ 


Integration of the simple integrals 




*~~ 


Recursion of type 


ml 




*~~ 


this is topMl 






*~~ 


Recursion of type 


m2 




*~~ 


this is topM2 






*~~ 


Recursion of type 


m3 




*~~ 


this is topM3 






*~~ 


Recursion of type 


m4 




*~~ 


this is topM4 






*~~ 


Recursion of type 


m5 




*~~ 


this is topM5 






*~~ 


perform integration 




*~~ 


Recursion of type 


tl 




*~~ 


this is treatnl 






*~~ 


Recursion of type 


nl 




*~~ 


this is treatnl 






*~~ 


Simplify 






*~~ 


Do the "rest"-inti 


3gration 




Time 


5.33 sec 


Generated terms = 


23 




d3179 


Terms in output = 


23 






Bytes used 


410 



d3179 = 

+ ep--3 * ( - 8/3 - l/3*a ) 

+ ep'-2 * ( 56/3 - 20/3*a ) 

+ ep--l * ( 112/3 - 16*z3 + 19/2*z2*a - 20*z2 - 97/12*a ) 

+ 334/3 + 1215/2*S2*a - 1620*S2 + 16*D3*a - 40*D3 - 1141/3*z3*a + 2368/ 
3*z3 + 144*z4*a - 288*z4 + 57*z2*a - 156*z2 - 32*a*B4 - 77/6*a + 64* 
B4; 

save problems/'PRB'/results/'DIAGRAM' .res 'DIAGRAM'; 
. end 
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